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ABSTRACT 

We present a detailed study of how the Star Formation Rate (SFR) relates to the interstellar medium (ISM) of 
M31 at ~ 140 pc scales. The SFR is calculated using the far-ultraviolet and 24 fim emission, corrected for the 
old stellar population in M3 1 . We find a global value for the SFR of 0.25 +o. 04 M© yr _1 and compare this with 
the SFR found using the total far-infrared (FIR) luminosity. There is general agreement in regions where young 
stars dominate the dust heating. Atomic hydrogen (H I) and molecular gas (traced by carbon monoxide, CO) or 
the dust mass is used to trace the total gas in the ISM. We show t hat the global surfa ce densities of SFR and gas 
mass place M31 amongst a set of low-SFR galaxies in the plot of Kennicutt (1998b). The relationship between 
SFR and gas surface density is tested in six radial annuli across M3 1 , assuming a power law relationship with 
index, N. The star formation law using total gas traced by H I and CO gives a global index of N = 2.03 ± 0.04, 
with a significant variation with radius; the highest values are observed in the 10 kpc ring. We suggest that this 
slope is due to H I turning molecular at Ecas ^ 10 M pc -2 . When looking at H2 regions, we measure a higher 
mean SFR suggesting a better spatial correlation between H 2 and SF. We find N ^ 0.6 with consistent results 
throughout the disk - this is at the low end of values found in previous work and argues against a superlinear 
SF law on small scales. 

Subject headings: Stars: formation, ISM: general, Galaxies: evolution, Local Group, Galaxies: spiral, Galax- 
ies: star formation, Infrared: galaxies, Submillimeter: ISM, Ultraviolet: stars 



1. INTRODUCTION 

Relationships between star formation and the quantity of 
gas in the interstel lar medium (ISM) seem to vary greatly 
across the universe ( Kennicutt & Eva nsl2012|) . There is a tight 
correlation between Star Formation Rate (SFR) and gas mass 
in most cases, but the slope of the relationship does not ap- 
pear to be consistent. A specified surface density of gas can 
correspond to a range of SF R even in a single galaxy. 

Schmidt's original paper (Schmidtl ll959l) . which focussed 
on the Milky Way, derived a relationship between the vol- 
ume density of star formation and the volume density of gas 

(Psfr ^ PgLs) w * m a P ower l aw index N = 2. The first ex- 
tragalactic measurements of the Schmidt law were carried out 
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bv lSanduleakl (11969b and lHartwickl (119711) on the Small Mag- 
ellanic Cloud (SMC) and M3 1 respectively and found indices 
ofiVsMC = 1-84 ±0.14 and N M31 =3.50 ±0.12. 

Since then, similar studies have tended to relate surface 
density of star formation (Hsfr) to surface density of gas 
(Leas), which are what we actually observe. However, the 
index probed using surface densities should be equivalent to 
that for volume densities as long as we have a constant scale 
height. 

A later study of 16 nearby galaxies b y lBoissier et al.l ([2003) 
found N = 2. IWong & Blitzl (120021) studied 6 nearby spi- 
rals and est i mated N to be in the range 1.2-2.1. However 
iHever et al.l (12004ft calculated an index of - 3.3 for M33 
when considering total gas, but N ~ 1.4 when looking only 
at m olecular hydrogen A more recent work on the same ob- 
ject (I Verley et al .1 120 1 ) find a wide range of indices (1.0 < 
N < 2.6) depending on gas tracer and fitting m ethod. In 
the co mprehensive (and most often cited) work of Kennicutt 
(1998b), he estimated the power index for 90 nearby galaxies 
using total gas (molecular gas only for starbursts) and found 
N = 1.40 ±0.15. A big question is whether this slope only 
works when considering global measurements or is it a mani- 
festation of a relationship on smaller scales. 

One interpretation of the Kennicutt result is that star 
formation timescales are dictated by the free-fall time, 
SFR ~M/r ff ~ pip- 1 ' 2 ~ p*' 2 , since r ff oc p' 1 / 2 
(e.g. lElmegreenl 119941 : iKrumholz & Thompson! 120071: 
Narayanan et al. 2008). Other work suggests that the super- 

*Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with important 
participation from NASA. 
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linear slope is a result of variations in the fraction of dense 
gas between normal spiral galaxies and starbursts and that the 
star f ormation law is linear given constant dense gas fraction 
(e.g. lLada et aT1l20l2h . 

Recently acquir ed data from the Gala xy Evolution Ex- 
plorer (GALEX, iMartin & GALEX Teamll2005l) and The H I 
Nearby Galaxies Survey (THINGS) has allowed the star for- 
mation law to be p robed on sub-kpc (~ 750 pc) scales as in 
Bigie l et all ([2008). They suggest that star formation is more 
directly related to molecular rather than total gas and find that 
the molecular gas star formation law follows a relationship 
with index, N = 1 (a linear relationship), consistently lower 
than the values they find for total gas. 

M31, due to its proxi mity and size (D = 785 kpc 
(McConnachi e et all |2005|) . apparent angular size is 1900 
gives the opportunity to probe these relationships over an 
entire galaxy, at smaller physical scales (comparable to the 
size of a giant molecular cloud) than all other e xtragalac- 
tic sp irals. It is also roughly solar metallicity (Yi n et al.l 
2009) so is a good analogue to the Milky Way. This has 
been done previously (on lower fidelity data than here) by 
Tabatabaei & Berkhuijsen ( 2010), who found a similar super- 
linear relationship between surface den sities of star form ation 
from Ha and total gas to that found in lKennicuttl d!998b|) for 
whole galaxies. 

Star formation tracers, whether looking at unobscured or 
embedded star formation, invariably rely on the assumption 
that the emission used as a SF probe originates dire ctly, or as 
a result of heating, from young stars (Calzetti 2007). This is 
a reasonable assumption in galaxies that have recently under- 
gone a starburst, as massive young stars burn brightly and die 
young, with less massive stars living much longer and provid- 
ing a minimal contribution to the ultraviolet lu minosity. How- 
ever, M3 1 has not undergone a recent starburst (Davidge et al. 
I2012h so contributions from older population s can have a sig - 
nificant effect on star formation estimates (|Calzettil |2012|) . 
This should, in principle, be possible to mitigate using trac- 
ers of the general stellar population. 

In this paper, we present multi- wavelength data of M3 1 and 
measure the total unobscured and embedded star formation 
rates separately using far-ultraviolet (FUV) and 24 /im data 
respectively. We compare this with the total gas, found by 
combining maps of neutral atomic hydrogen (H I) and carbon 
monoxide (CO( J=l-0)), which traces the molecular hydrogen 
(H 2 ). 

The maps tracing SFR and gas mass are used to calibrate the 
SFR and gas mass found using the far infrared (FIR) emission 
from M31, as obs erved with the Herschel Space Observa- 
tory dPilbratt et al.ll2010|) as part of the Herschel E xploitation 
of Lo cal Galaxy Andromeda (HELGA) project (Frit z et all 
120121) . We compare our SFR from UV and 24 /am emis- 
sion with that found from FIR luminosity. The interstellar 
gas mass is also traced using the dust mass estimated from 
the FIR Spectral Energy Distribution (SED), scaled using the 
observed gas to dust ratio. Here we aim to see how well this 
gas map correlates with SFR, hence whether dust mass traces 
star forming regions. 

Finally, we use this collection of SFR and gas maps to probe 
the power law relationship between SFR surface density and 
the gas surface density, or Kennicutt- Schmidt (K-S) law. Our 
analysis is performed on individual pixels in M3 1 and investi- 
gates how the law varies with different gas tracers on sub-kpc 
scales. 



2. DATA 

Our first method of tracing star formation uses 
GALEX FUV and NUV observations of M31 
(IMartin & GALEX TeamI 120051) . alon g with warm dust 
emission seen in Spitzer MIPS 24 /im ([Gordon et al. 2006) 
and s tellar emission from Spitzer IRAC 3.6 /im ( Barmbv ^et al.l 
2006) (FigureCQ. 

The HELGA collaborati on obtained obser vations of M3 1 
in five Herschel bands (jFritz et al.l 12012b . They are 
PACS jPoditschet al.1 HoiQh 100 and 160 /im and SPIRE 
dGriffinetal] 12010b 250. 350 and 500 /im. Details of the 
data reductio n for both P ACS and SPIRE maps can be 
found in [F ritz et al] d20l2h . The Spitzer MIPS 70 /im map 
dGordon et al.ll2006|) . is employed to extend the wavelength 
range for our calculation of the FIR spectral energy distribu- 
tion (Figure [2]). 

We independent ly probe the interstellar medium using H I 
dBraun et al.ll2009h and CO( J=l-0) maps dNieten et al.f 2006) 
(Figure 0. Note that the CO map covers a smaller area than 
the H 1. The area not covered by the CO will be used in the 
calculation of total gas, providing there is sufficient H I. 

In this work, we divide the maps into elliptical annuli of 
constant galactocentric radius. We do this to test the effect of 
radius on the star formation law, with the option to relate this 
to the Toomre Q criterion, which relates to rotational velocity 
and shear. It also allows us to isolate the 10 kpc ring, where 
the majority of star formation in M31 is occurring; and the 
central regions which are dominated by an older stellar pop- 
ulation. The right-hand image in Figure [3] shows the colour 
coding of datapoints used in all subsequent plots depending 
on their radial distance from the centre. The ellipses are cre- 
ated assuming a po sition angle of Andromeda of 38° and an 
inclination of 77° dMcConnachie et all 12005). The colours 
indicate how the datapoints from those regions will be repre- 
sented in later figures. Distances are in units of Rm m which 
we take to be 21.55 kpc (Ide Vaucouleurs et aDl991[) . 

For analysis, the maps are individually smoothed and re- 
gridded to three pixel scales, based on the lowest resolution 
map used in the analysis. We modify the full- width half- 
maximum (FWHM) beamwidth to match the effective Point 
Spread Function (PSF) by Gaussian smoothing the image us- 
ing the IRAF function imgauss. The maps are regridded using 
the IDL astrolib function, FREBIN. Any offsets in the coordi- 
nates of the pixels are corrected using wcsmap and geotran in 
IRAF. 

The first scale used here is the highest resolution star forma- 
tion map we can create using the FUV and 24 /im emission 
as a tracer. This corresponds to the lowest resolution (MIPS 
24 /im) FWHM beamwidth of 6" (a heam = 2.55 ") and a 
pixel size of 1.5 " . This scale is applied to the 3.6 /im, 24 /im, 
NUV and FUV maps. 

We aim to study the relationship between SFR and gas mass 
on the smallest scales attainable. To this end, we also use 
maps smoothed to the resolution and grid size of the neutral 
atomic hydrogen map, again the lowest resolution map used 
here. The effective FWHM beamwidth is 30 " (abeam = 
12.7") with a 10 " pixel size. This scale is applied to the data 
mentioned above, with the addition of the CO( J=l-0) map. 

In order to compare gas mass (from H I and CO( J=l-0)) 
and star formation in M3 1 (from FUV and 24 /im emission) 
with the Herschel observations, the majority of the analysis 
is performed on a scale corresponding to the beamsize of 
the lowest resolution SPIRE map (500 /im). These images 
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Fig. l . — : Images used in the creatio n of the FUV and 24 /im star formation ma p of M31. From left: GALEX FUV and N UV 
maps (Martin & GALEX Teamll2005l) . Spitzer IRAC 3.6 /im (iBarmbv et al.ll2006l) . SpitzerMIPS 24 /xm dGordon et al.ll2006l) . 
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Fig. 2.—: Top, from l eft: Spitzer MIPS 70 /im dGordon et al.ll2006l) . PACS 100 and 160 /im dFritz et al.ll20l2h : bottom: SPIRE 
250, 350 and 500 /im (iFritz et al.H20l2h . The PACS and SPIRE observations are from the HELGA collaboration. 



have an effective FWHM beamwidth and grid size of 36 " kernels from lBendo et al.l (1201 2|) as described in Smit h et al. 
(o"beam = 15.5 "). Since the beamwidth and pixel size are (l2012al) . 
equivalent, the pixels can be described as approximately 'in- 
dependent,' as there is no correlation between them. Here, 
the MIPS maps (Figure [2]) were smoothed using convolution 
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Fig. 3.—: Left, integrated H I emission dBraun et al.l l2009): centre, CO( J=l-0) dNieten et al.ll2006l) : right, colour key of elliptical 
annuli. Each annulus is 0.2 i?M3i thick apart from the inner two which are 0.1 i?M3i thick. 



3. STAR FORMATION RATE 

3.1. FUVand24/j J m 

The star formation rate is first calculated from the GALEX 
F UV and Spitzer 24 a m maps, using the method prescribed 
in iLeroy et aD (|2008|) . However, to expand on this we also 
use GALEX NUV and Spitzer IRAC 3.6 fim maps to correct 
for foreground stars and emission from old stellar populations 
respectively. 

FUV emission is predominantly from unobscured high- 
mass stars (O, B and A-type), so this tracer is sensitive to star 
formation on a time s cale of ~ 100 Myr (e.g.|Kennicutt 1998a; 
ICalzetti et all 120051: ISalim et alJ 12007b , 24 /mi emission is 
predominantly due to dust-heating by UV photons from bright 
young stars, and is sensitive to a star formation timescale 
of < 10 Myr (e.g. [ Calzetti etaD 120051: iPerez-Gonzalez et alJ 
2006; ICalzetti et alJl2007h , 

The star for mation surface den sity is calculated using the 
formulation in ILeroy et ail (120081) which uses a Chabrier ini- 
tial mass function (IMF): 



A second problem is that some of the emission could be 
from an older stellar populat ion. This is a genera l problem and 
not specific to M31 (e.g. iKennicutt et aiTl2009h . We expect 
this to be a bigger iss ue near the centre of the galaxy. P revious 
FIR work on M31 (Tabatabaei & Berkhuijsen 2010) avoids 
this problem by measuring the SFR at radii greater than 6 kpc 
only, based on the assumption that the centre of the galaxy 
contains negligible star formation. Old stars are fainter but 
redder, so emit relatively stronger at 3.6 fim. This means we 
can mitigate for the old stars by determining 2fuv /h.6 (here- 
after, olfuv) an d h^/h-Q (hereafter, 0^4) in regions where 
we assume star formation has ceased, and use this to remove 
the component of FUV and 24 fim emission coming from old 
stars. So, the emission we associate with star formation is 
given by, 

^FUV,SF = ^FUV — <^FUV ^3.6 (2) 



'24,SF 



a 2 A I[ 



3.6- 



(3) 



Ssfr = 8.1 x 10- 2 7 FUV + 3.2±J-2 x lO" 3 1 



(1) 



where Hsfr has units of M yr -1 kpc -2 and FUV and 
24 fim intensity (/) are in MJysr -1 . The pixel size corre- 
sponds to a distance of ~ 140 pc. If comparing like for like 
with other galaxies, an inclination correction factor of cos i 
(where the inclination of M31, i = 77°) must be included 
in order to 'deproject' the image, effectively giving values 
as they would be for a face-on galaxy. This prescription as- 
sumes all the 24 fim emission in M3 1 is due to dust heating by 
newly formed stars, and that the FUV is emitted exclusively 
by young stars. There are, of course, other sources of these 
tracers which are unrelated to star formation which must be 
taken into account. 

The first issue is foreground stars. These are se l ected and 
removed using the UV colour, as in ILeroy et aD (120081) — 
if Inuv/^fuv > 15 the pixel is blanked in both the FUV 
and 24 am map ( some 24 fim emission will be stellar, e.g. 
Bendo et al. 2006). We assume this ratio will only be reached 
where a pixel is dominated by a single star, which, given our 
resolution will never be associated with M3 1 . 



ILeroy et all (12008b explored this by looking at the ratio of 
fluxes determined in elliptical galaxies. They found apuv = 
3 x 10 ~ 3 and 0^24 = 0.1. However, if we compare the 
3.6 fim emission with the FUV and 24 fim in M31, we see 
that these values are not necessarily appropriate here (Figure 
[4]). The 24 fim emission in the bulge (shown by red points) 
follows the ratio found in ellipticals (black-dashed line, Fig- 
ure |U right), so we will use the same value for 0^4. c^fuv 
is found to be much lower here (Figure [4] top). We speculate 
this is due to dust extinction in M31, which is not an issue 
in passive elliptical galaxies as they con tain little dust (e.g. 
Smith et a l. 2012b: Rowlands et al.ll2012h . It is also stated in 



Leroy et al.l (120081) that there is a large scatter in this ratio so 
a discrepancy is not surprising. An independent correction is 
found by performing linear fits on the inner regions of M3 1 . 
Ellipses within a radius 0.05, 0.1 and 0.2 i?M3i give gradients 
(a F uv) of 8.42 x 10~ 4 , 7.99 x 10~ 4 and 7.44 x 10~ 4 respec- 
tively. Here, the mean value, apuv = 8-0 x 10 -4 , will be 
employed to correct FUV emission for the old stellar popula- 
tion in M3 1 . We performed this analysis on the high resolu- 
tion maps to maximise the number of datapoints and checked 
that the slope was consistent with that found using the lowest 
resolution (36 ") maps. 
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Fig. 4.—: In order to study the effect of old stars on our star formation tracers, we plot FUV vs 3.6 /im and 24 /im vs 3.6 /im 
emission, in this case on the high resolution map. Colours indicate the distance, from the galactic centre, of each p ixel (see 
Figure 0. The black dashed trendline indicates the correction for the old stellar population used in iLeroy et al.l (|2008|) . based 
on /fuv/^3.6 and I24/I3.6 found in ellipticals. The solid trendline is the best fit to FUV vs 3.6 /im (left) and 24 /im vs 3.6 /im 
(right) in the inner regions of M31 (r < 0.1i?M3i)- We can see that there is a tight correlation in both plots indicating that FUV 
and 24 /im emission in the centre of M31 is predominantly from old stars. 



Once this correction is applied, we have a map of surface 
density of star formation, Esfr in units of M yr _1 kpc -2 
(Figure[5]). The correction in this work has the effect of re- 
ducing the measured global SFR by ~ 25%. 

It should be noted that when looking at the region imme- 
diately outside the 10 kpc ring only (dark blue points, behind 
red points (not visible), Figure |U left), we see a tight cor- 
relation between the FUV and 3.6 /im emission. Unlike the 
centre, however, there are pixels in that region that do not fol- 
low this correlation. This indicates that despite a significant 
population of old stars in the ring, star formation is still oc- 
curring at significant rates compared to the rest of the galaxy 
(-O^Moyr- 1 ). 

We can scale the star formation map to give the to- 
tal star formation occurring within a pixel. This allows 
us to determine the global SFR for M31, which we find 
to be 0.25±g;g2 

M yr x , almost an order of magnitude 
lo wer than the canon i cal v alue of ~ 2M yr _1 , quoted 
in IChomiuk & Povichl (1201 lh for the Milky Way but con- 
sistent with the lower limit of ~ 0.27 M yr -1 found in 
iTabatabaei & Berkhuijsenl (120101) for M31. 

3.2. Star formation from far-infrared luminosity 

Star formation can also be calculated using the total FIR lu- 
minosity. This ideally probes the embedded SFR and is sensi- 
tive to cooler dust temperatures. This can be an issue in deter- 
mining total SFR in dust-deficient galaxies where significant 
starlight is not attenuated by dust but should not be a prob- 
lem here. An issue that is relevant to M31 is that a significant 
component of the d ust heating could be from an evolved stel- 
lar population (e.g. iBendo et al.1 1201(1 iBoquien et al.l 1201 it 
iBendo et alJbOia Smit h et al.ll2012al) . 

The total FIR flux was integrated in frequency space using 
a linear interpolation between the six datapoints (70-500 /im) 
for each pixel independently. Each value was converted from 
a flux to a luminosity in L assu ming a distance to M31 of 
785 kpc (iMcConnachie et al.ll2005l) . 

If we assume FIR luminosity is exclusively re-radiated 



emission from warm dust that is heated during a continuous 
starburst, the FIR luminosity is equal to the total luminosity 
of the starburst. The total SFR, Msf is then, 



M SF = S M F (W/lO iU L ) M yr 



(4) 



where Smf depends on the assumed IMF of the region be- 
ing studied and the timescale of the starburst. Changing 
these assumptions gives radically different conversion factors. 
For example, assuming a Kroupa IMF with a star formation 
timescale of 10 Gyr gives 5mf ~ 0.6 wher eas the same IM F 
with a timescale of 2 Myr gives S M f ~ 3.2 dCa lzetti 20 O)- 

Here we employ the value quoted in iKennicu tt (1998b|) of 
^mf = 1-7, which assumes a Salpeter IMF with a low mass 
cut-off of 0.1 M and a timescale of ~ lOOMyr. We should 
state here that this assumes a continuous starburst which keeps 
consistency with the previous method. This conversion factor 
gives a global star formation rate of 0.52 M yr _1 . 

If, as before, the old stellar population has a significant 
effect on the dust heating at these wavelengths, we would 
naively expect to see a correlation between the FIR luminosity 
and 3.6 /im emission. However, the total FIR luminosity is a 
function of dust mass and dust temperature, so if the distribu- 
tion of dust is different t o the distribution of stars (as it is in 
M31. [Smith et al.ll2012a|) there will be no correlation, even i f 
the old stars are the major heating source (Ben do et~aTll2012[) . 
This was tested and no correlation is visible. 

Without any kind of correction for the old stars, the star for- 
mation rate from the FIR emission is measured to be approx- 
imately doub le the estimate from the FUV and 24 /im tracers 
(Section [37Tb . This is expected as M31 has not gone through 
a starburst in its recent history, so a significant portion of the 
heating is due to the interstellar radiation field (ISRF). 

As discussed in the previous section, past work on M31 
elected to omit the central region of the galaxy when deter- 
mining the global SFR, due to the dominance of old stars in 
this region. If we omit the central region out to 0.2 Rm31 the 
measured SFR reduces from 0.52 M yr _1 to 0.48 M yr _1 . 
This minimal difference suggests that the over-estimate is not 
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limited to heating from old stars in the bulge. This is consis- 
tent with the correlation obser ved b etween 3.6 /im emission 
and 24 fim in the ring (Section[3j]), indicating old stars have 
a significant heating effect here also. 

3.3. Comparison of star formation tracers 

The star formation maps made using FUV and 24 fim emis- 
sion, and that from FIR luminosity can be seen in Figure [5] 
The FUV and 24 /im tracer has units of M yr _1 kpc -2 , but 
we have elected to display the FIR star formation tracer in 
terms of FIR luminosity as the conversion factor between lu- 
minosity and star formation rate is very uncertain. 

Sub-mm wavelengths are more greatly affected by heating 
due to the ISRF than the 24 fim emission, as this regime is 
sensitive to cooler dust temperatures. This means the FIR 
emission is susceptible to heating from more distant stars 
(those that are not in the same pixel) , ma king determination of 
a correction factor difficult (Section |3^2| >. In order to compare 
the two SFR tracers, we modify the FUV and 24 /im tracer 
to match the assumptions made in creating the map of Ssfr 
fr om FIR lum i nosity . 

iLeroy et all (|2008[) state th at their tracer retur ns a star for- 
mation rate 1.59 x lower than lKennicuttl (Il998bl) . from where 
we get our 5my (Equation [4]). Because of this, in Figure[6]we 
com pare Ssfr from FIR luminosity (as described in Section 
13.2b with the same from FUV and 24 /im with an additional 
factor of 1.59 to approximate the same IMF, and no correc- 
tion for old stars. The low-mass cut-off (mi = 0.1 M ) of 
the IMF and star formation timescale (r = 100 Myr) are al- 
ready equivalent. This gives a global SFR ofO.51M yr _1 , 
consistent with the value from FIR luminosity. 

It does appear that the SFR as measured from FIR lumi- 
nosity is slightly lower relative to the FUV and 24 /im tracer 
in the very centre and outer regions of M31 (Figure [6] right). 
There is a possible issue with PACS observations not recov- 
ering all of the flu x in low surface brightness regions (e.g. 
Ania no etall 12012). If this was the case, it may contribute 
to the observed effect between annuli but our global SFR will 
be minimally affected. We do not believe this discrepancy is a 
major issue here as this map is not used in any further resolved 
analyses of low surface brightness regions. 

Despite the general consistency between SFR tracers, we 
should remain aware that since the conversion factor between 
tracer luminosity and star formation rate depends on the IMF, 
we are not necessarily recovering the correct value. It is pos- 
sible that the IMF we assumed for M3 1 is not appropriate, or, 
because the star formation rate is low, we are not sampling 
the whole IMF leading to fluctuations in the tracer luminosity 
for a fixed SFR. This can be an issue for a variety of tra cers, 
including IR and UV emission ( Kennicutt & Evans 20 l3). 

For the analysis that follows, we elected to use the com- 
bined FUV and 24 /im emission as our star formation tracer, 
as we are able to correct for the old stellar population. We ar- 
gue that this gleans more reliable SFRs in low SFR regimes, 
like those in M3 1 , than when using the FIR luminosity. 

4. THE INTERSTELLAR MEDIUM IN M3 1 

The ISM is made of predominantly neutral atomic and 
molecular hydrogen. A map of total gas can be produced by 
summing these two constituents and multiplying by a factor of 
1.36 to account for heavier element abundances (mostly He- 
lium), or alternatively b y assuming the t otal gas is well traced 
by dust emission (e.g. lEales et al .1120121) . 



4. 1 . Total gas from H I and CO observations 

The Hi map is taken from iBraun et al.l (|2009). In order 
to keep consistency with our maps of star formation, and to 
allow comparison of galaxies with different inclinations, i, we 
employ a factor of cos i to 'deproject' the galaxy. 

H2 is the most abundant molecule in the ISM, but lacks 
a dipole moment so is not easily observable. For this rea- 
son, CO (usually J=l-0) is empl oyed as a tracer mo lecule. 
For M31 we use the map from iNieten et ail (120061) . The 
conversion between CO emission and quantity of H2 is still 
a contentious topic and uses the so-called X -factor (e.g . 
| Wallll2007b iGloyer & Mac Lo^l201ll: iNaravanan et al.H 201 U 
Feldman n et al.H201 lblBolatto et al.H2013l) . where~ 

7V H2 /cm- 2 = X x Ico/Kkmr 1 (5) 

The conversion factor specific to molecular clouds in 
the north eastern arm of M 31 was argued to be 5.68 x 
10 2Q (K kms" 1 !" 1 cm" 2 in ISofue etall (11994 . This was 
found by estimating virial masses (from their size and veloc- 
ity width) and comparing to the CO line intensit y. This is 
larger than the value found by iBolatto et al.l (12008) of ~ 4 x 
10 20 (Kkrns -1 )" 1 cm -2 using the same method. However, 
the ISM of M3 1 is dominated by neutral atomic hydrogen so 
it is not clear whether the virial masses provide an overesti- 
mate of the mass of molecular hydrogen in these clouds. Here 
we will assume X = 2 x 10~ 2Q (K kms" 1 )" 1 cm" 2 (e.g. 
Strong et al. 1988; Pineda et al. 2010), which agrees with the 
value derived in ISmith et al.l (l2012al) . Any constant discrep- 
ancy in the X -factor will result in a horizontal translation in 
our logio(SH 2 ) vs logio(^SFR) plots and so will have no ef- 
fect on the calculation of our K-S index for molecular gas. It 
may skew the calculation using total gas but the effect is likely 
to be small due to the dominance of H I in the ISM of M3 1 . 

At this point we should note the su ggestion that metallicity 
has a n effect on the X -fact or (e.g. Ilsraell [T997I: IStrong etall 
12004 . ISmith et al.l (I2012al) found a radial variation in the 
gas-to-dust ratio suggesting a metallicity gradient in M3 1 and 
hence a gradient in X. Any variation should not be a big issue 
when looking at total gas, as the ISM in M31 is dominated by 
neutral atomic hydrogen. We do need to keep this uncertainty 
in mind when dealing with molecular gas, however. 

4.2. Total gas traced by dust 

The interstellar gas can, in pri nciple, also be tra ced by the 
distribution of dust in a galaxy (Eales et al. 2012). The dust 
map of M31 is taken from Smi th et al.l (120 12a!) / where dust 
mass is found by fitting a modified blackbody function in each 
pixel where there is a 5 a detection in all bands. This should 
mitigate a gain st the low-surface brightness issues discussed 
in S ection 13. 3 1 

In lSmith et all d20l2al) a fit was performed to gas-to-dust vs 
radius, to determine how the conversion factor varies. It was 
found that the relationship is logarithmic with radius varying 
between ~ 30 near the centre and ~ 100 in the 10 kpc ring, 
consistent with the value found in the MW (Spitzer 1978). We 
use this function to create a second total gas map as traced by 
dust. In the following section, as with the other ISM tracers, 
this will be used to observe how well dust mass correlates 
with star formation. 
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Fig. 5. — : Left, star formation rate map obtained using the FUV and 24 jim emission (FWHM beam width = 6" ', pixel size = 
1.5 "); right, total FIR luminosity at A > 70 fim (FWHM beamwidth and pixel size = 36 "), found by integrating the emission 
from Herschel and Spitzer observations shown in Figured The SFR from infrared luminosity is found using a constant conversion 
factor based on the assumed IMF and assumed length of the continuous starburst. 
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IO gio( £ SFR(,R, /M so,y r " 1 ) r/R M31 

Fig. 6.—: Left: Esfr (star formation rate surface density) found from FIR luminosity (assuming a Salpeter IMF) vs Hsfr from 
FUV and 24 fim emission (scaled to match values assuming a Salpeter IMF and without a correction for the old stars). The 
colours represent the radius of the pixel from the galactic centre (Figure [3}. The solid black line indicates a 1:1 relationship, 
the dashed black lines indicate factors of 4 offsets. Right: ratio of Ssfr from FIR luminosity to Ssfr from FUV+24 fim star 
formation surface density with radius. The errorbars represent one standard deviation of the scatter in this value across each 
elliptical annulus. 
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5. THE STAR FORMATION LAW 

In this section, we probe the star formation or K-S law as- 
suming the following relationship, 

£ S FR = -4£(L;, ( 6 ) 

where TV is the power index and A is related to the star for- 
mation efficiency (SFE). 

We separately look at how M31 compares to other local 
galaxies in terms of global SFR surface density (calculated 
from FUV and 24 /im emission) and gas surface density; and 
what relationship the star formation law follows on a pixel- 
by-pixel basis when considering various components of inter- 
stellar material. 

5.1. Global star formation law 

Figure [7] compares the mean surface density of star 
formation rate with the mean surface density of gas for 
globa l measurements of galaxies from Kennicutt (1998b) 
and lLeroy et al.l (|2008), with corresponding glob al values for 
M31 overplotted. The SFRs from this paper and Lerov et al. 
( 2008) are scaled to match the assumptions made in Kennicutt 
( 1998b). The mean values for M31 are found over all pixels 
with sufficient signal-to-noise in both maps (SFR and gas). 
The difference in measured SFR is due to the different selec- 
tion effects depending on the gas tracer. We can immediately 
see that the mean SFR for regions containing sufficient H2 
(ico > 5 aco) or dust (/ > 5 a in five Herschel bands) is 
higher, suggesting a better spatial correlation between SFR 
and both molecular hydrogen and dust tha n total gas. 

The low gas mass galaxies studied by iKennic utt (1998b) 
generally appear to have higher star formation rates than M3 1 , 
although we note that they estimate star formation using a 
different SF tracer. However, early-type spirals like M31 
are expected to ex hibit a low SFR per unit area, as stated in 
iKennicuttl d!998a|) . Mean surface densities of both total and 
molecular gas are consistent with the same parameters for nor- 
mal spirals studied in previous work. 

5.2. Resolved star formation law 

From our gas and SFR maps, we should be able to inves- 
tigate the Kennicutt-Schmidt star formation law, on a 'per 
pixel' basis across the galaxy. This section aims to test how 
calculation of this law changes with different gas tracers. It 
has been suggested that H 2 is a bette r tracer of star formation 
than total gas (e.g. Bigi el et aT1l2011|) . although it is not clear 
if this is the case in M31 where H I dominates the ISM. 

Here we use data from our maps of surface density of star 
formation as found from FUV and 24 /im emission, against 
surface density of total gas, molecular hydrogen only and gas 
traced by dust. We select pixels that satisfy Ssfr > 5 ctsfr 
and Seas > 5 acas, where ctsfr is the standard deviation 
of the background of the star formation map and acas is a 
combination of the uncertainties of the constituent gas maps 
(e.g. for total gas this will be the scaled uncertainties in the 
integrated H I and CO( J=l-0) images). 

5.2.1. Fitting 

Here we perform a linear fit in order to find the index, N 
from equation [51 assuming that 

logio S S fr oc TV logio S Gas . (7) 

In Figure [51 the signal-to-noise (S/N) cuts are clearly man- 
ifest. In total gas, the major cut-off is horizontal (limited 



by ctsfr); in molecular gas the cut is vertical (limited by 
ctho)- Previous works appear to exhibit a similar cut-off (e.g. 
Tabatabaei & Berkhuijsen 2010) but with no attempt to mit- 
igate for this when performing a fit. After some exploration 
we conclude that the signal to noise cut does indeed bias the 
data and must be mitigated against (see Appendix). We test 
other methods of fitting and find that binning the data in or- 
der of increasing star formation gives the most reliable return 
gradient when testing the relationship to total gas. 

We therefore attempt to mitigate for the S/N cut in our data 
using the same method. When looking at the total gas from H I 
and CO measurements, we order in bins of increasing SFR, 
with an equal number of datapoints (500) in each. We then 
plot the mean surface density of gas (I^Gas / M kpc _2 )in 
each bin, against the mean surface density of SFR (Ssfr / 
M kpc _2 yr _1 ) and perform the fit on these points in the 
logarithmic domain using a least squares routine in MATLAB. 

In the case of H2 only, the S/N cut-off is more apparent in 
gas mass so we bin the data in order of increasing gas mass, 
with 100 points in each bin. 

Gas mass estimated from dust mass exhibits a more com- 
plex selection effect so binning is not attempted here. The 
majority of the points that are omitted correspond to points 
that appear towards the low- SFR regime of the bottom-left 
window in Figure [U This will affect our calculation of the SF 
law but we believe the analysis is still valid as we are pref- 
erentially selecting regions that are more important for star 
formation. 

5.2.2. Kennicutt-Schmidt index 

Figure [8] shows plots of SFR versus gas mass (in units of 
M pc -2 , to keep consistency with previous work) for each 
pixel, with trendlines for each annulus (top row) and a global 
fit (bottom row). Figure [9] shows K-S index as a function of 
radius for the various gas tracers. Global values are indicated 
by a dashed line and are also given in Table Q] 

The K-S index for each annulus varies between 1.0 and 2.3 
when considering total gas, with the higher values applying 
to the lOkpc ring. The global value for total gas is ~ 2.0, 
decreasing to ~ 1.8 if we double the X-factor. 

The star formation law with H2 gives a shallower gradient 
(N = 0.6) but is more constant between annuli. Doing the 
same with gas traced by dust gives a similarly shallow slope. 

5.3. Discussion 

Our star formation law in M31 using total gas gives a K-S 
index N ~ 2.0, significan t ly hig her than the value found by 
Tabatabaei & Berkhuijsen (2010) but consistent with values 
found in previous work on other galaxies. One possible expla- 
nation for a steep slope is H I becoming optically thick. How- 
ever, when performing the analysis on an opacity corrected 
map of H I, we see no real change in the calculated star for- 
mation law. An alternative is that the hydrogen turns molec- 



TABLE 1 

Global K-S parameters using different tracers. 



Gas tracer 


N 


<j n 


log(A) 


a a 1 dex 


R 


Total gas (Hi + CO) 


2.03 


0.04 


-19.11 


0.31 


0.6 


H2 only from CO 


0.60 


0.01 


-7.41 


0.08 


0.7 


Total gas from dust 


0.55 


0.01 


-7.55 


0.07 


0.7 



Note. — Ssfr = A SQ as . The quoted standard deviation refers to that found using 
the fit on binned data. The correlation coefficient, R refers to raw data points. 
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O K98, normal spiral, total gas 

O K98, normal spiral, H2 only 

* K98, starbursts, H2 only 

x L08, HI only 

V L08, H2 {upper limit only) 

x L08, H2 only 

x L08, Total gas 

This work, M31 , total gas 
^ This work, M31, H2 only 

This work, M31 , gas from dust 




log /M DC 

3 10 x gas sol r 



Fig. 7.—: Global SFR vs gas mass, derived in M31 using three gas tracers, with Kennicutt (1998b) and iLeroy et alJ (l2008|) 
galaxies. G a laxies are plotted using total gas, H I only and H2 only, dependin g on availability of data. SFRs from this work and 
ILeroy et ail (|2008|) are scaled to match the assumptions of Kennicutt (1998b). The das hed diagonal lines represent the gradient 
the galaxies should follow given a Schmidt law of the type found in Kennicutt (1998b) (N = 1.4, black) or a linear relationship 
(N = 1, blue), as has been suggested by several recent papers (e.g. Rahman et al. 20 lob- 



ular (which occurs at Seas ~ lOM pc -2 ) but is not traced 
by our CO( J=l-0) map in these high SFR regimes, hence t he 
turnover apparent in previous work (e.g. Bigie l et aDl2Q08|) is 
not visible here. 

We get a lower mean star formation rate when looking at all 
regions of gas which suggests we are not isolating star form- 
ing regions as well as when using other gas tracers. The po- 
sitions of datapoints for each pixel in Figure [8] (bottom row) 
for total gas, appear to depend on their radius. This is consis- 
tent with the fits to each annulus showing a horizontal offset 
(top row), with the inner regions to the left of the plot (low 
gas mass). This would suggest that the threshold for star for- 
mation changes with radius (see also Figure 0). However we 
should note that the inner regions contain relatively fewer dat- 
apoints. Also, it is possible that despite our correction, we still 
overestimate star formation in the centre. 

The index found here using molecular gas (TV ~ 0.6) does 
argue against a superlinear relationship on small scales. It is 
possible the resolution of our images may affect calculation 
of our index but when testing the relationship using a 500 pc 
pixel scale, we see a negligible difference. Some recent work 



has suggested that the star formation la w is linear (N= 1) 
when looking at molecular clouds (e.g. Rahman et al. 2012|) 
and that the superlinear relationship of Kennicutt (19983) is 
not a manifestation of a relationship that applies on smaller 
scales, but may be the result of systematic differences be- 
tween the galaxies that are related, but not limited to, gas mass 
alone. 

A 'sub-linear' relationship (N < 1), indicates that star for- 
mation is less efficient at high gas densities which would be an 
intriguing result. We urge caution with our value for the index 
with molecular gas however, due to the significant scatter. 

One other issue we should keep in mind is that because the 
scales we are probing are small compared to other extragalac- 
tic sources, the surface density of star formation and surface 
density of gas may not be directly relatable to the correspond- 
ing volume densities as our scale height is more likely to vary 
between regions. 

6. SUMMARY 

In this paper we have determined the surface density of star 
formation in M31 using combined FUV and 24 jim emis- 
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Fig. 8.—: Star formation rate surface density against gas surface density for three gas tracers. The top row shows fits to the 
Kennicutt-Schmidt law for each annulus. The bottom row shows the data for every pixel studied, with a global fit indicated by a 
solid black trendline, with the calculated index N quoted in the top left corner of each panel along with the correlation coefficient 
of unbinned data, R. Each column denotes a single gas tracer. From left, they are: total gas from H I and CO( J=l-0), assuming 
a CO-Ho conversion factor of 2 x 10~ 20 (Kkms -1 )- 1 cm , with an additional factor of 1 .36 for heavier element abundances: 
molecular hydrogen (H2) from CO( J=l-0); and total gas traced by dust mass (see Smit h et aLll20 12a). assuming a radial gradient 
in the gas-to-dust ratio. 



sion and separately the far-infrared luminosity. We aim to 
correct the former for emission from both unobscured and 
embedded old stars and find a global star formation rate of 



0.25^0.04 M© yr _1 . The FIR emission appears to be corre- 
lated with the SFR map made using FUV and 24 /im emis- 
sion. However, we are unable to correct for the old stellar 
population as there is no correlation visible between FIR lu- 
minosity and 3.6 fim emission in the galactic centre. 

We produce two maps of the total gas in M31. The first 
uses Hi and CO, assuming a CO-H2 conversion factor of 
2 x 10 2Q (Kkms -1 )" 1 cm" 2 . We use the r adially varying 
gas-to-dust ratio found in lSmith et al.l (1201 2a|) to produce the 

second map of total gas from the dust emi ssion. 

When compar i ng wi th previous work by Kennicutt (1998b) 
and iLeroy et al.l (|2008) on the global SFR and gas mass, we 
find the mean molecular gas surface density and star forma- 
tion rate surface density fo r M31 sit on the low end of the 
relation determined in lKennicuttl (Il998b|) . 

Our measurement of the star formation law on sub-kpc 



scales varies with gas tracer. The most direct measurement, 
using H I and CO to trace total gas, gives power law index 
N ~ 2.0 when looking at the whole galaxy, consistent with 
the range of values found in previous work but we believe this 
slope is a result of H I saturation. The values measured in ra- 
dial annuli vary between 1.0 and 2.3, with the highest values 
being measured in the 10 kpc ring, where the vast majority of 
star formation is occurring. 

Using molecular gas only gives a much lower K-S index 
of N ~ 0.6, suggesting that a superlinear relationship with 
molecular gas is not applicable on sub-kpc scales in M3 1 . 
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APPENDIX 
FITTING 

Our UGas against Hsfr plots exhibit clear signal-to-noise (S/N) cut offs which it is imperative we address. The first task was 
to test whether the cuts result in a bias when performing the fit. 

We created two arrays, x and y where x contains all integers between —1000 and +1000 and y = mx where — 1 < m < +3. 
Gaussian noise is applied to both the x and y values to simulate the observed spread in points. We then apply a cut at a specified 
y value, again mimicking the data (Figure [TOb. 

The poly fit algorithm in MATLAB is used to perform the fit on the data above the cut. Poly fit is a least- squares routine that 
minimises residuals in the y-axis parameter. When m ^ the calculated gradient is consistently shallower than the input, 
indicating a bias. We attempt to mitigate for this by ordering the data in bins of increasing y with an equal number of points in 
each bin. We replace this data with a single point based on the the mean or median of the binned data. The fit is then performed 
using the same algorithm on these averaged points. In both cases we get slightly steeper gradients when — 1 < m < +1 
converging to near perfect agreement in the region where 1.5 < m < 2.0. In all cases studied, the unmitigated fit is more deviant. 

We therefore attempt to mitigate for the S/N cut in our data using the same method. When looking at the total gas from Hi and 
CO measurements, we order in bins of increasing SFR, with an equal number of datapoints (500) in each. We then plot the mean 
gas mass in each bin, against the mean SFR and perform the fit on these points, using the logarithmic units. 

In the case of H2 only, the SNR cut-off is more apparent in gas mass so we bin the data in order of increasing gas mass, with 
100 points in each bin. Gas mass as estimated from dust mass exhibits a more complex selection effect so binning is not attempted 
here. 
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Fig. 10.—: Simulated Kennicutt- Schmidt index fitting. Left: Example simulated dataset with input gradient of 1.0. Black points 
are the selected data, cyan points are those that have been discarded before performing the fit. Right: Input gradient versus 
measured gradient for a range of input gradients and fitting methods. In both plots, green represents the unmitigated fit; blue 
represents the mean of binned data; red, the median. The black trendlines represent the input gradient. 



